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ABSTRACT 

This paper extends our previous study of planet /brown dwarf accretion by giant 
stars to solar mass stars located on the red giant branch. The model assumes that the 
planet is dissipated at the bottom of the convective envelope of the giant star. The 
giant's evolution is then followed in detail. We analyze the effects of different accretion 
rates and different initial conditions. The computations indicate that the accretion 
process is accompanied by a substantial expansion of the star, and in the case of 
high accretion rates, hot bottom burning can be activated. The possible observational 
signatures that accompany the engulfing of a planet are also extensively investigated. 
They include : the ejection of a shell and a subsequent phase of IR emission, an increase 
m the 7 Li surface abundance and a potential stellar metallicity enrichment, spin-up of 
the star due to the deposition of orbital angular momentum, the possible generation 
of magnetic fields and a related X-ray activity due to the development of shear at the 
base of the convective envelope, and the effects on the morphology of the horizontal 
branch in globular clusters. We propose that the IR excess and high Li abundance 
observed in 4-8% of the G and K giants originate from the accretion of a giant planet, 
a brown dwarf or a very low-mass star. 

Key words: accretion - planetary systems - stars:low-mass, brown dwarfs - stars:latc 
type - stars:evolution - starsrmass loss - stars: chemically peculiar - stars:rotation - X- 
ray:stars - infrared:stars 



1 INTRODUCTION 

Planets and brown dwarfs have now been detected around 
several main sequence stars (e.g. Mayor & Queloz 1995, But- 
ler & Marcy 1996, Cochran et al. 1997, Rebolo et al. 1995, 
1996, Basri et al. 1996, Marcy et al. 1998, Butler et al. 1998, 
Delfosse et al. 1998), and strikingly, a fraction of these "hot 
Jupiters" orbit very close to the central star, at less than 1 
AU. This proximity to the central star poses a theoretical 
problem regarding the mechanism that brought these plan- 
ets to these orbits, but it also raises the question of their 
fate. Indeed, Rasio et al. (1996) showed that the orbit of 51 
Peg's companion is tidally unstable and that the planet will 
ultimately be engulfed in the stellar envelope. The problem 
we address here is that of the effects of the accretion of a gi- 
ant planet or a brown dwarf on the structure and evolution 
of solar mass stars on the red giant branch. This study is 
the continuation of our previous work (Siess & Livio 1998), 
in which we investigated the swallowing of massive plan- 
ets/brown dwarfs by a 3M@ AGB star. 

The general scenario for the accretion process is as fol- 



lows : planets in close orbits (as observed) will be engulfed 
in the envelope of giant stars as the latter evolve away from 
the main sequence. Due to viscous and tidal forces, angu- 
lar momentum is imparted to the envelope and the planet 
spirals-in. Previous studies of the spiralling-in process (e.g. 
Livio & Soker 1984; Soker, Harpaz & Livio 1984) have re- 
vealed that depending on the mass of the planet/brown 
dwarf three different evolutionary paths can take place. If 
the mass of the planet is larger than some critical value 
Merit, the planet can accrete a fraction of the envelope mass. 
On the other hand, lower mass planets will either evapo- 
rate or collide with the stellar core. For the specific stellar 
model used in their computation, Livio and Soker (1984) 
estimated M cr i t ~ 20Mj up . However, these results must be 
taken with great caution since the physical processes were 
treated only approximatively. We also emphasize that such 
computations are very sensitive to the stellar model, no- 
tably to the size and mass of the convective envelope. More 
recently, Sandquist et al. (1998) performed 3D hydrodynam- 
ical simulations of Jupiter-like planets impacting onto solar 
type main sequence stars. They also found that depending 
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on the initial stellar model, planets can be evaporated or 
survive the crossing of the convective zone. 

The present work has been motivated by the fact that 
planets surrounding low mass stars are more likely to be en- 
gulfed in the star's envelope during the red giant branch 
(RGB) phase than during the asymptotic giant branch 
(AGB) phase. Indeed, if we compare (for different stars) 
the ratio of the maximum radius a star can reach during 
the RGB to its maximum radius on the AGB phase, this 
ratio is much closer to one for low mass stars (M ^ 2Mq) 
than more massive stars (e.g. Soker 1998a). Consequently, 
we expect that in most cases low mass stars will interact 
with close low-mass binary companions on the RGB. 

In the present paper, we used a spherically symmetric 
stellar evolution code to follow the evolution of a solar type 
star that accretes a giant planet or a brown dwarf on the red 
giant branch. In we describe briefly the computational 
approach and the initial models. Then, in we analyze 
the effects of accretion on the structure and the evolution of 
the star. We present results for different mass accretion rates 
and different initial stellar models on the RGB. The observa- 
tional signatures expected from this process are investigated 
in detail in §0. A discussion and conclusions follow. 



2 INITIAL CONDITIONS 

To investigate the effects of planet accretion by solar type 
stars, we used 3 different evolutionary models of the sun cor- 
responding to 3 different points on the RGB. These models 
are characterized by their mass, which is equal to M = 1.0, 
0.9 and O.8M and they are referred to as model #A, #B 
and #C, respectively. The main properties of the initial 
models are summarized in Table |l[ The general RGB struc- 
ture is a degenerated He core, surrounded by a thin hydro- 
gen burning shell (HBS) that provides most of the luminos- 
ity. A thin radiative layer separates the HBS from the large 
convective envelope. The star was assumed to lose mass as 
prescribed by Reimer's law (1975) at a rate given by 

Moss = -3.98 x 10~ 13 n R -^ Mq yr" 1 , (1) 

where the parameter r/R — 0.5 in our computations, and L, 
R, and M are given in solar units. 

The input physics for the accretion process has been 
previously described by Siess and Livio (1998). Briefly, to 
determine the locus of the dissipation process we estimated 
the Virial temperature of the brown dwarf and the critical 
radius at which tidal effects lead to a total disruption of the 
brown dwarf. For characteristic values of the mass and radius 
for the brown dwarf, we found a temperature of a few million 
K and a critical radius of the order of ~ 0.1 — IRq. From the 
numbers presented in Table [j], we can localize the accretion 
process at the bottom of the convective envelope near the 
HBS. To estimate the accretion rate M acc , we determined 
the orbital decay timescale of the brown dwarf in the vicinity 
of the dissipation region. This timescale roughly corresponds 
to the time required by the brown dwarf to cross the dis- 
sipation region and thus represents the rate at which mass 
is added there. For typical stellar parameters, we inferred 
mass accretion rates of the order of 10 -4 — 10~ 5 Mq yr _1 , 
depending mainly on the mass of the brown dwarf. 



Table 1. Physical properties of the initial models 





#A 


#B 


#c 


Mass (Mq) 


0.9928 


0.9004 


0.8072 


Radius (-Rq) 


10.72 


88.71 


137.10 


L (L ) 


45.2 


1133.6 


2016.3 


age t (10 9 yr) 


12.1398325 


12.1999969 


12.2022346 


fc t nv 


0.3324 


0.2366 


0.2098 


M enVibot (Mq) 


0.2583 


0.3904 


0.4336 


M H BS,top (Mq) 


0.2474 


0.3883 


0.4322 


M corc (Mq) 


0.2419 


0.3867 


0.4310 


Rcnv.bot (Rq) 


0.543 


0.791 


0.831 


R-HBS.top (-R©) 


0.08586 


0.1033 


0.1060 


Rcorc (Rq) 


0.02587 


0.02538 


0.02492 


Pcnv.bot (gem -3 ) 


4.057 10" 2 


1.586 10" 3 


8.715 lO" 4 


/"HBS.top (gem" 3 ) 


2.733 


0.3726 


0.2510 


Pcorc (gem -3 ) 


369.7 


135.2 


112.0 


Tenv,bot (10 6 K) 


3.286 


2.324 


2.209 


THBS.top (10 6 K) 


11.0 


12.3 


12.6 


T corc (10 6 K) 


29.8 


44.5 


48.2 


T c (10 6 K)* 


34.6 


60.9 


69.9 



tgyration radius [(7 cnv /M cnv i? 2 nv )l/2] 
* central temperature 



During the accretion phase, we adopted a uniform mass 
deposition profile as described in Siess and Livio (1998). 
More specifically, we deposit in each numerical shell an 
amount of accreted matter proportional to the mass of the 
shell, starting from the top of the HBS and up to a radius 
r = IRq. Our simulations indicate that small changes in 
the profile of mass deposition do not substantially affect the 
results. The main point here is the depth at which the mass 
is deposited inside the star, and in the following sections 
we identify the accretion zone as the region located close 
to the base of the convective envelope where most of the 
energy accompanying the accretion process is released. Fi- 
nally, if the convective envelope reaches the HBS (beginning 
of Hot Bottom Burning or HBB), the matter is deposited 
from the base of the convective envelope (which is inside 
the HBS). In our computations, we also assumed that the 
matter is deposited with the same entropy as the one of 
the local matter, implying that the rotational energy has 
been used to heat up the accreted matter (e.g. Harpaz & 
Soker 1994). The deposition of accreted matter is also ac- 
companied by modifications of the chemical abundances. We 
adopted a solar metallicity with relative abundances given 
by Anders and Grevesse (1989) for the accreted matter and 
we assumed that deuterium had been previously burnt in the 
brown dwarf. Finally, for numerical convenience, we turned 
off the mass loss during the accretion phase and used grey 
atmosphere models to speed up the convergence process. 



3 NUMERICAL RESULTS 

The stellar evolutionary code used in these computations 
has been described in detail in several papers (e.g. Fores- 
tini & Charbonnel 1997, Siess et al. 1997) and we refer 
the reader to those publications for more details. Compu- 
tations are presented for 2 characteristic accretion rates of 
10 -4 and 10 _5 MQyr _1 , and they involve different masses 
for the planet/brown dwarf (see Table 0). 
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Table 2. 


Numerical models 






Models 


A^ini (M Q ) M acc 


(Moyr" 1 ) 


accreted mass 


AO 


0.99 


10 -5 


0.10 


Al 


0.99 


10~ 5 


0.01 


A 2 


0.99 


10~ 4 


0.10 


A3 


0.99 


10~ 4 


0.01 


BO 


0.90 


IO" 5 


0.10 


131 


0.90 


IO" 5 


0.01 


B2 


0.90 


IO" 4 


0.0051 


B3+ 


0.90 


io~ 4 


0.0058 


CO 


0.81 


IO" 5 


0.10 


CI 


0.81 


IO" 5 


0.01 


C2 


0.81 


IO" 4 


0.0010 


C3 


0.81 


10- 4 


0.0025 


C4+ 


0.81 


IO" 4 


0.0030 



^For these models, we were not able to follow reliably the evolu- 
tion and the computations were stopped at the indicated accreted 
mass 



1000 




4 6 

log (t -t ) (yr) 

Figure 1. The evolution of the radius Ft and luminosity L 
for the different models computed with a lower accretion rate 
(Afacc = 10 -5 Mq yr — 1 ). The time ip corresponds to the age of 
the initial model and is indicated for the different models in Table 
[l The solid lines refer to the standard solar model (no accretion), 
the dotted and dashed lines correspond to models in which the 
accreted mass is equal to 0.01 and O.IMq, respectively. The de- 
picted models are AO, Al, B0, Bl, CO and CI. 



3.1 Evolution of the structure for lower accretion 
rates 

The deposition of mass in the star is accompanied by the 
release of a large amount of gravitational energy that acts 
to expand the envelope (Fig. ^). However, because of the 
long timescale of heat diffusion, the expansion is delayed, 
especially for the low luminosity stars (cases AO and Al). 
As the radius increases, the star also becomes more lumi- 
nous. According to Reimer's law, these modifications to the 
structure significantly increase the mass loss rate, and we 
expect some changes in t he s ubse quent evolution and stel- 
lar environment (see fjhjand §4.5). The expansion proceeds 



until a thermal equilibrium is reached, in which the rate of 
energy loss by the surface balances the rate of energy de- 
position by the accretion process. Depending on the initial 
model, the return to contraction occurs after ~ 6200, ~ 470 
and ~ 200 yr, in cases AO, B0 and CO, respectively. This 
timescale corresponds roughly to the time of heat diffusion 
(thermal timescale) of the expanding part of the star and it is 
approximatively given by r t horm ~ GMM env /2RL lux , where 
M onv is the mass of the envelope and Lace is the accretion 
luminosity. The latter is given by Laccfr) = GM T Maxx/r, 
where A/ acc is the mass accretion rate and M r and r the 
mass and radial coordinates of the accretion zone. It should 
be noted that L acc represents an upper limit to the accre- 
tion luminosity since energy is lost during the spiralling-in 
process to heating and spinning up the envelope. However, 
throughout most of the envelope, the orbital energy disspa- 
tion rate represents a small fraction of the intrinsic stellar 
luminosity and thus we do not expect that the stellar struc- 
ture will be significantly affected until the planet reaches 
the accretion zone. Finally, we should mention that the en- 
ergy deposition might not be local (in the accretion region) . 
For example, the star could be tidally excited, oscillations 
may be initiated and then damped elsewhere in the enve- 
lope. Also, in this common envelope phase, magn etic fields 
might be generated in the accretion region (see £4.3) and 



energy could be dissipated at the surface of the star or in 
the corona (e.g. Regos & Tout 1995). 

During the envelope expansion, the opacity k increases 
in the accretion region and rapidly the radiative gradient 
becomes larger than the adiabatic one (V ra d oc kL > V a d); 
thus the convective envelope deepens (Fig. ^ . Thereafter, as 
the HBS advances in mass and reaches the cool regions that 
have been previously expanded, the nuclear energy produc- 
tion rate decreases significantly. As an order of magnitude, 
during most of the accretion phase, the luminosity due to 
H burning represents only a few percent of the total stel- 
lar luminosity. The accretion luminosity thus supports the 
convective envelope almost entirely and the accretion region 
acts as a source shell that prevents an efficient H burning. 

To see how the structure is affected by the accretion 
process, we have presented in Fig. ^ the profiles of several 
variables for case AO, but they are globally similar in the 
other models. In the accretion region where a large amount 
of energy is deposited, the luminosity and temperature rise. 
But, due to the relatively long timescale associated with heat 
diffusion, energy accumulates locally and both profiles now 
present a bump. The heat excess finally diffuses when the 
star has sufficiently expanded and a monotonically increas- 
ing luminosity profile is restored. However, the temperature 
profile still presents a local maximum that persists through- 
out the accretion phase. This is due to the fact that this 
region undergoes a compressional heating from the mass ac- 
cumulated above it and also from the expanding region sur- 
rounding the core. Indeed, the luminosity L r is negative just 
behind the luminosity jump, indicating a region in expan- 
sion (e grav < 0). The te mpe rature is thus compelled to rise 
and, as we shall see in §3.3, this situation can give rise to a 
thermal instability. The temperature gradient also becomes 
progressively steeper beyond the bump (Fig. []) as energy 
accumulates. This contributes to the shrinkage of the HBS 
and to the subsequent drop in the nuclear luminosity. 

In the HBS, prior to accretion, the nuclear energy 
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Figure 2. The evolution of the boundaries of the HBS and the 
convective envelope. The dashed lines delineate the burning shells 
(enuc > 5 ergg - 1 s _1 ) and the dotted line indicates the locus of 
maximum energy generation. The hatched areas indicate convec- 
tive regions. The different models are labeled in each panel and 
were computed with an accretion rate equal to 10 — 5 M© yr — 1 . 
Note the extinction of the HBS in case AO (lower right panel). 
In the upper right panel (case CO), a thermal instability devel- 
ops and leads to the development of a convective tongue. The 
dotted-dashed lines delineate the He burning shell. 
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Figure 3. The evolution of the profiles of luminosity (L r ) and 
temperature (T) for case AO. The solid line refers to the initial 
model, the dotted, short-dashed and long-dashed lines correspond 
to models for which 0.019, 0.039 and 0.086AfQ have been accreted, 
respectively. The accretion rate is equal to 10 — 5 A/q yr — . The 
peak in the lower left panel represents the nuclear energy produc- 
tion rate in the initial model and thus delineates the H burning 
shell. 



production rate is mainly supplied by the CNO bi-cycle. 
More precisely, the most important reactions (and their 
contribution) are 13 C(p,7) 14 N (~ 25%), 14 N(p, 7 ) 15 (~ 
23%), 15 N( P ,q) 12 C (~ 16%) and ls O(/3+, v) 15 N (~ 16%). 
However, because the HBS advances in mass into a cooler 
region, the ppl branch progressively takes over the nuclear 
energy production and during the accretion phase, the reac- 
tions 2 3 He(7,2p)a and 13 C(p,7) 14 N contribute ~ 45% and 
~ 10% of the total nuclear energy production, respectively. 
On the other hand, in the region of maximum temperature, 
13 C is now efficiently burnt through the 13 C(«,n) 16 reac- 
tion. This reaction is highly energetic and provides up to 
~ 20% of the global nuclear energy production (which still 
remains very low, however). In Fig. W, we show the evolu- 
tion of the nuclear activity for case B0. First, 13 C burning 
is ignited in the vicinity of the temperature peak (dotted 
line). Then, as it is depleted and converted into le O, the nu- 
clear energy production rate progressively decreases near the 
temperature peak (short-dashed and long-dashed curves), 
and finally we end up with a configuration in which 13 C is 
burnt through 13 C(a,n) 16 around the temperature peak 
and through 13 C(p,7) 14 N above the bump in the HBS. In 
case CO, the temperature in the peak reaches ~ 1.2 10 8 K 
and He is ignited off-center in a non d ege nerate shell (we 
will discuss this special case further in §3.3). 



3.2 Higher accretion rates 

For higher accretion rates, the situation is somewhat dif- 
ferent because of the appearance of hot bottom burning 
(HBB). Indeed, contrary to the lower accretion rate case, 
the larger energy release contributes to a further increase in 
the opacity and the convective envelope deepens even more. 
It finally reaches the HBS where its advance stops (Fig. |j). 
Our computations indicate that HBB is activated approxi- 
mately 200, 25 and 15 yr after the onset of accretion in cases 
A2, B2 and C2 respectively. The temperature at the base of 
the convective envelope rises suddenly from ~ 2 — 3 to about 
13— 18 million Kelvin, depending on the initial model. How- 
ever, this temperature is still too low to give rise to a rich 
nucleosynthesis. The onset of HBB entails a significant and 
fast expansion of the envelope as heat is directly transferred 
to the convective zone. Moreover, the arrival of 3 He from the 
convective envelope and its destruction through the reaction 
2 3 He(7,2p)a increases temporarily the nuclear luminosity in 
the HBS. During the HBB phase, 7 Li is slowly depleted be- 
cause of the competing effects of Li production through the 
reactions a( 3 He,7) 7 Be(e~, ^) 7 Li and Li destruction through 
the reaction 7 Li(p, a)a (see also §4.1). 

Finally, in cases B2 and C2 we were not able to follow re- 
liably the computations. The star's tendency to expand and 
increase its luminosity made the computations intractable. 
At that point the envelope is probably ejected and the star 
enters a proto-planetary nebula phase. 



3.3 The relaxation of the star 

After the end of the accretion phase, the envelope rapidly 
contracts because of the sudden suppression of accretion en- 
ergy. However, in the models in which the duration of the 
accretion process is shorter than the thermal timescale as- 
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Figure 4. Evolution of the nuclear energy production rate (e n uc), 
temperature (T), and 16 and 13 C chemical abundances for 
model BO (M a cc = 10 — 5 Mq yr -1 ). The solid lines refer to the 
initial model t = to. The dotted, short-dashed and long-dashed 
lines correspond to an age t — to equal to 6791, 9760, and 11040 
yr, respectively. Note the persistence of the temperature excess 
during the accretion phase. 



sociated with the envelope (as in cases Bl and CI for exam- 
ple) , the expansion still proceeds for a while after the end of 
mass addition. Thereafter, the gravitational energy released 
during the contraction phase heats up the central parts of 
the star and the HBS re-ignites. During that period of con- 
traction, the convective envelope retreats towards the sur- 
face and HBB disappears. The contraction rate is very fast, 
especially for stars having a shallower convective envelope. 
The computations indicate that in a few thousand years, the 
radius decreases by a factor of ~ 2 — 3 which, we suspect, 
could initiate stellar pulsations. When the HBS finally pro- 
vides most of the stellar luminosity (<; 80%), the contraction 
ends and the star resumes a standard evolution, climbing the 
red giant branch until He flashes at its tip. Interestingly, the 
star is now more luminous and has a larger radius than in 
the standard evolution despite the small amount of added 
mass (Fig. |l|). The enhanced mass loss rate that follows can 
affect the subsequent evolution (see §4.5). 

In all our simulations, except for case CO, the tem- 
perature "bump" dissipated as the gravitational energy re- 
leased during the contraction phase heated up the stel- 
lar interior. However in this particular model, in which 
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Figure 5. The evolution of the internal structure (left and lower 
right panels) and radius (upper right panel) for the high accretion 
rate case (A-f acc = 10 — 4 Mq yr - 1 ). The evolution of the structure 
is depicted in a similar way to that in Fig. ^. Note the appearance 
of hot bottom burning in these computations. In the upper right 
panel, the solid, dotted, short-, long- and dot-dashed lines refer to 
the standard case and to cases A2, B2, C2 and C4, respectively. 



the accretion process is maintained for ~ 10 4 yr at a rate 
equal to 10 -5 Mq yr -1 , the temperature in the peak reaches 
~ 1.2 10 8 K and He ignites off-center in a non degenerate 
shell. The ignition of He through the 3a reactions is very 
energetic and this gives rise to a thermal instability (Fig. 
[j) . The geometrical thinness and the large opacity of the He 
burning shell (HeBS) are responsible for the thermal run- 
away that leads to the development of a convective tongue. 
The release of energy due to He burning is used to expand 
the layers above and secondarily to heat up the core. The 
HBS is pushed outward into a cooler region and it extin- 
guishes. The remaining energy excess is then absorbed by 
the envelope which subsequently expands. During the ther- 
mal pulse, the nuclear energy production comes from 3a 
(~ 70%), 14 N(a,7) 18 F (~ 17%) and, to a lesser extent, from 
ls O(a, 7 ) 22 Ne and ls O(a, 7 ) 22 N (~ 7% each). The thermal 
pulse is also responsible for the further destruction of H 
through 12 C(p,7) 13 N. After a significant expansion that al- 
lows the star to get rid of the energy excess, contraction 
resumes. However, due to the extremely low H content left 
after the pulse (the H mass fraction X < lO -20 ), the HBS 
cannot ignite before a new thermal instability is triggered. A 
new He shell flash thus occurs deeper inside the star where 
both the temperature and He abundance are higher. As evo- 
lution proceeds, the pulse duration increases but the max- 
imum luminosity due to He burning decreases mainly be- 
cause of the progressive lifting of the core degeneracy as the 
temperature increases. During the subsequent evolution, the 
star experiences relatively short phases of expansion follow- 
ing the thermal pulses and longer phases of contraction. The 
deepening of the pulse accounts for the smaller and smaller 
increase in radius following the He shell flashes since most of 
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log (t-g (yr) 

Figure 6. The evolution of the radius, He luminosity and internal 
structure for case CO. In the upper panel, the different lines follow 
the same convention as in Fie 



the released energy is absorbed in the radiative layers before 
reaching the convective envelope. After the 12 th pulse, the 
HBS ignites precisely at the mass coordinate correspond- 
ing to the top of the first convective tongue in which the 
H abundance is high (X 3> 10 -20 ). The H luminosity then 
increases and the thermal pulse finally reaches the center of 
the star. Because of the previous thermal instabilities, the 
degeneracy has been lifted in the interior and consequently 
He does not flash. Helium is now burnt at the center and the 
star will spend the next ~ 4 10 7 yr on the Horizontal branch. 



4 POTENTIAL OBSERVATIONAL 

SIGNATURES OF AN ACCRETION EVENT 

In this section, we attempt to give a broad overview of the 
different observational signatures that might be expected 
from the accretion of a planet or a brown dwarf by a giant 
star on the red giant branch. In general, effects may be ex- 
pected in at least six different areas, which we describe in 
some detail below. The modifications to the structure dur- 
ing the accretion process increase the mass loss rate and may 
lead to the ejection of shells. They also have potential effects 
on the morphology of the horizontal branch. The deposition 
of angular momentum during the spiralling-in process can 
also alter the stellar structure by increasing the rotation 
rate and this may also produce X-ray emission. Finally, the 
mixing of the dissipated planet material in the envelope can 
change the surface chemical abundances. All of these effects 
may have important observational consequences but, as we 
shall see in the following sections, it is more the combination 
of several signatures, such as a high rotation rate and/or a 
large IR excess and/or a high lithium abundance that can 
identify an accretion event less ambiguously, rather than the 
detection of one particular signature. 



4.1 Modifications of the photospheric chemical 
composition 

During the red giant phase, chemical changes are mainly the 
result of "dilution" effects due to the deepening of the con- 
vective envelope but light elements such as 7 Li, 9 Be and X1 B 
can also be burnt as the convective envelope approaches its 
deepest extent. The first dredge up reduces the 7 Li abun- 
dance by more than 2 orders of magnitude. There is also 
a decrease in the 12 C/ 13 C ratio by a factor of 3, changing 
from about 83 to 26. The first dredge up also accounts for 
the large enrichment of 3 He as a by product of H burning. 

With our assumed chemical composition for giant plan- 
ets and brown dwarfs (from Anders & Grevesse 1989), the 
deposition of accreted matter, mostly in the convective en- 
velope, will principally modify the chemical abundances of 
the light elements that were previously depleted during the 
first dredge up. In that respect, 7 Li will be the most affected. 
Our computations indicate that, depending on the accreted 
mass, the 7 Li abundance can be increased by more than 2 or- 
ders of magnitude, making this e lem ent easily detectable at 
the surface of the staij^ (see also §4.6). As an illustration, we 



present in Fig. [?] the 'Li surface abundances resulting from 
our computations and compare them with what is expected 
from simple dilution effects. In the dilution hypothesis, we 
assume that the envelope mass M onv is held constant during 
the accretion process and that no nuclear burning is taking 
place in the convective zone. The mass fraction Xi of an 
element i is then simply given by 



Xi 



xr v x Monv + xr x Ma, 

Mace + M onv 



(2) 



where M aC c is the mass accreted into the envelope, and Xf nv 
and Xf° c are the mass fraction of element i in the enve- 
lope and in the accreted matter, respectively. Note that this 
expression is general in the sense that Xi can either rep- 
resent the H mass fraction X, the He mass fraction Y or 
the metallicity Z. Figure ^ shows that the simulations fol- 
low relatively well the theoretical curves defined by Eq. (H) . 
The slight deviations from this relation are mainly due to 
nuclear burning that causes Li depletion. Note that in case 
A2, 7 Li is severely depleted due to the long lasting HBB. 
The deposition of accreted matter will also increase the 9 Be 
and 11 B chemical abundances, but to a lower extent, since 
the effects of dilution during the first dredge up are smaller 
for these elements. A modest increase in the 12 C/ 13 C ratio 
can also take place. 

The effects of planet accretion on the stellar metallic- 
ity are very sensitive to the amount and chemical compo- 
sition of the accreted material and also to the mass of the 
convective envelope (which depends on the evolutionary sta- 
tus of the star). If the metal content of the giant planet is 
higher than the one adopted here, the metallicity can be 
further enhanced. Secondly, if the planet is engulfed while 
the star is still on the main sequence for example, the con- 
vective envelope is then very shallow (Af on v — O.O2M0 for 
the sun), and consequently the stellar photospheric compo- 
sition can be substantially modified (Eq. Pj). Recent deter- 



* Following de la Reza et al. (1996), we define as a "Li rich G or 
K giant" a star with a Li abundance larger than log e (Li) = 1.2 
[where by definition loge(H) = 12.00] 
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Figure 7. The evolution of the 7 Li surface abundance as a func- 
tion of the accreted mass. The solid lines represent the 7 Li abun- 
dance as a result of dilution for the 3 initial envelope masses 
characterizing our initial models (Af cnv = 0.74, 0.51 and 0.37AfQ 
in models A, B and C, respectively). The dotted lines represent 
the results of our computations for cases AO, B0, CO and A2. Note 
the large Li depletion in case A2 due to the long lasting HBB. 



minations of the chemical composition of the main sequence 
stars hosting planets (Gonzales 1997, 1998) indeed show that 
in some cases the metallicity is abnormally high, suggesting 
that proto-planetary material or planets have been accreted 
(e.g. Lin et al. 1996, Laughlin & Adams 1997). More re- 
cently, Sandquist et al. (1998) showed that a Jupiter-like 
planet can be evaporated in the convective envelope of the 
sun. They also inferred that in some cases a metal enrich- 
ment as high as 50% can be obtained, but the exact value 
depends greatly on the mass of the convective zone and on 
the chemical composition and interior structure of the giant 
planet. 

All the above suggests that if the presence of planets 
around solar type stars is rather common, and if some of 
these planets are indeed in close orbits, then certain chemical 
peculiarities in giants stars may be the consequence of the 
accretion of planets. 



4.2 Rotation 

During the spiralling-in process, angular momentum is im- 
parted to the envelope. As a consequence, the stellar rotation 
rate will be substantially increased. More quantitatively, if 
all the orbital angular momentum of the planet is entirely 
deposited and redistributed in the envelope, i.e. if we as- 
sume that the convective envelope is rotating rigidly, then 
the ratio of the envelope angular velocity u to the critical 
velocity (surface Keplerian angular velocity) uik is given by 

£-Hfe)(3?r(£r(fr- <»> 



where a is the orbital separation of the brown dwarf, k the 
gyration radius and M env the mass of the envelope. Equa- 
tion (^) indicates that only massive planets or brown dwarfs 
(Afbd > 5Mj up ) can spin- up the giant envelope to a signifi- 
cant fraction of the Keplerian velocity (J; 10%). The spin-up 
would be most effective at the tip of the giant branch, where 
both Menv and k are smaller. 

Measurements of the projected rotational velocity of lu- 
minous class III giants indicate a sudden drop in vsin i near 
spectral type GO III (Gray 1989). The same discontinuity 
is also observed in subgiants at spectral type F8 IV (Gray 
& Nagar 1985) and on the main sequence near mid-F stars 
(Kraft 1970). These observations thus delineate in the HR 
diagram a rotational dividing line (RDL) where slow rota- 
tors (v sin i ^ 5 kms -1 ) populate the cool side (de Medeiros 
et al. 1996b). The strong braking at the rotational disconti- 
nuity occurs very rapidly, perhaps on a time scale as short 
as ~ 10 5 yr or less (Ayres et al. 1998) , and probably results 
from an efficient magnetic braking attributed to a dynamo 
action (Gray 1981, 1982). In the "rotostat" hypothesis (Gray 
1986), angular momentum is dredged up from the rapidly 
rotating core by the deepening of the convective envelope. 
This activates the dynamo and then an efficient magnetic 
braking takes place. In this model, the rotation rate is con- 
trolled in a self-regulatory way by the balance between the 
arrival of angular momentum from the core and the strong 
braking that ensues, and consequently, the rotational ve- 
locity is maintained constant. However, a few stars do not 
conform to this scheme and still rotate rapidly beyond the 
RDL, challenging the rotostat hypothesis. For example the 
case of 7 Boo = HR 5225 remains unexplained (Gray 1989, 
1991). This single star of spectral type G5III has a moderate 
lithium abundance log£(Li)=1.2 and its rotational velocity 
(vsin i = 14.1 kms -1 ) is three times higher than the average. 
This isolated case could fit into a scenario in which the giant 
star has recently accreted a planet and its surface layers have 
not yet spun down. The moderate lithium abundance would 
suggest the accretion of a relatively small planet. More gen- 
erally, de Medeiros et al. (1996a) found that among a sample 
of ~ 900 stars, less than ~ 5% of the late G and K giants 
located to the right of the RDL have v sin i > 5 kms - . 

On the other side (left) of the RDL, FK Comae-type 
stars show extreme rotation velocities, (w/wk <; 0.10). The 
general explanation for their fast rotation assumes that these 
single G or K giants are undergoing (or have just completed) 
coalescence with a close binary companion (e.g. Webbink 
1977, Bopp & Stencel 1981, Livio & Soker 1988, Welty and 
Ramsey 1994). Other scenarios have also been proposed 
and they generally invoke tidal synchronization in a very 
close binary, normal evolution from a rapidly rotating up- 
per main sequence star, or the dredge up of angular mo- 
mentum (e.g. Simon & Drake 1989). But, as pointed out 
by Rucinski (1990), FK Coma has a few times more (> 8) 
angular momentum than a typical rapidly rotating main se- 
quence star. Therefore, some angular momentum must have 
been added (accreted) to the star. Generally FK Comae-like 
stars lack photospheric lithium detection. However in some 
cases, e.g. HD 36705 (Rucinski 1985), HD 33798 (Fekel & 
Marschall 1991) or IE 1751+7046 (Ambruster et al. 1997), 
a high lithium abundance has been observed (see also §p~(|). 
These observations clearly support our accretion scenario 
since they combine two different outcomes of the accretion 
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process : a high rotation rate and a high Li abundance. How- 
ever, different scenarii have also been proposed and among 
them the possibility that these stars were initially fast ro- 
tating A or F type stars that have rapidly evolved in the 
Hertzsprung gap. The suggestion is then that they recently 
became giants and are now depleting lithium due to the 
deepening of the convective envelope (Fekel 1988). It has 
also been suggested (Fekel & Marschall 1991) that angular 
momentum and CNO processed material can be dredged up 
from a rapidly rotating core as the convective zone deepens. 
In this scenario Li can be produced and burnt through the 
chain 3 He(a, 7) 7 Be(e, v) 7 Li(p, a) 4 He and, if on its way to the 
surface Li is not entirely depleted, then we may account for 
temporarily high rotation together with enhanced lithium 
abundance. However, this hypothesis would imply that es- 
sentially every star arriving to the RGB should go through a 
phase of fast rotation and high Li abundance. The paucity of 
such candidates poses a severe problem to this scenario and 
can be regarded as additional evidence against the rotostat 
hypothesis. 

To conclude this section, we feel that observations sup- 
port the coalescence hypothesis for the explanation of FK 
Comae-like stars and we stress that rapidly rotating stars 
may also have accreted angular momentum from the swal- 
lowing of a giant planet or a brown dwarf (see also Viten- 
broek et al. 1998). The planet accretion scenario can also 
account for enhanced Li abundances and in that respect, 
HD 33798 and HD 36705 are good candidates for such an 
event. Finally, we would like to point out that fast rotation 
velocities, of the order of 15-40 kms -1 , have also been ob- 
served in horizontal branch stars (e.g. Peterson et al. 1995, 
Cohen & McCarthy 1997). Notably, Peterson et al. (1983) 
already mentioned that the deposition, during the red gi- 
ant phase, of the orbital angular momentum of a planet of 
a few Jupiter masses could account for the fast rotational 
velocities observed in horizontal branch stars. Some AGB 
stars also exhibit high rotation rates (e.g. V Hydrae). For 
these exceptional stars, the rapid rotation rate is thought to 
be due to spinning up by a companion that has been swal- 
lowed or is in a common envelope phase (Barnbaum et al. 
1995). More generally, the deposition of angular momentum 
in the envelope of evolved stars by a companion (even sub- 
stellar) is frequently invoked to explain the axisymmetrical 
morphology of bipolar nebulae (e.g. Iben & Livio 1993, Livio 
1997). 

4.3 X-Ray emission 

If the accretion of angular momentum can cause the star to 
rotate differentially then differential rotation and convection 
can generate a magnetic dynamo. The magnetic field can 
then confine the plasma in the stellar chromosphere and 
produce X-ray emission activity. To see if this scenario is 
plausible, let us first examine the observational situation. 

Observations indicate the existence of an X-ray dividing 
line in the HR diagram (XDL also called coronal dividing 
line) which delimits giant stars with (to the left) and with- 
out coronal emission characteristic of high temperature gas 
(T > 10 K). In the HR diagram, this line is approximatively 
vertical at spectral type ~K3 and it spans the luminosity 
classes II, III and IV (Haisch et al. 1991, 1992). Note that 
bright giants of luminosity class II and lb turn out to be 



X-ray sources if exposed sufficiently deeply (Reimers et al. 
1996). On the right side of the XDL, the vast majority of the 
stars exhibit a very low X-ray emission (at least below the 
actual detection threshold) whereas a large spread of X-ray 
luminosities is observed for G and K giants of earlier spectral 
type than K3 (Strassmeier et al. 1994, Hiinsch et al. 1996). 
The X-ray emission originates from the stellar corona and 
requires a strong magnetic field in order to confine the hot 
plasma, since the thermal velocity of the particles exceeds 
the escape velocity. Changes in the magnetic field configura- 
tion (due to changes in the dynamo) are generally invoked to 
account for the sudden drop in X-ray emission at the XDL 
(e.g. Rosner et al. 1991, 1995) but the precise reasons why 
the dynamo is modified remain unclear. Recently Hiinsch 
and Schroder (1996) presented a "revised" XDL that seems 
to follow the evolutionary track of a 1.2Mq star. If this is 
confirmed, then evolutionary effects may account for this 
demarcation line. Low mass stars, on the right of the XDL, 
are generally slow rotators on the main sequence and conse- 
quently they sustain a lower magnetic activity. Conversely, 
more massive stars (M J; I.6M0) have faster rotation rates 
on the main sequence and therefore they can develop an ef- 
ficient dynamo and produce a higher X-ray emission. More 
recently, Hiinsch et al. (1998) looked for X-ray emission in 
the "forbidden" region beyond the XDL and they detected 
4 intrinsically X-ray bright M-type giants : 15 Tri, HR 5512, 
42 Her and HR 7547. No evidence for binarity exists yet for 
15 Tri and HR 5512 but even if the presence of a compan- 
ion is confirmed, it cannot easily explain the X-ray emission 
which remains quite unusual^. In a ROSAT survey of a com- 
plete, volume limited sample of late type stars, Hiinsch et 
al. (1996) also pointed out the unusual behavior of S Sgr = 
HR 6859, an apparently single, slowly rotating but strong 
X-ray emitting K3 III giant. Another puzzling case is that 
of HR 1362 (Strassmeier et al. 1990). Despite its low rota- 
tion rate this single G8 IV giant presents an X-ray emission 
that is more than one order of magnitude higher than that 
expected from the empirical rotation versus activity trend. 
The strong emission of HR 1362 led Stepien (1993) to sug- 
gest that it might represent an evolved Ap star that has 
preserved its strong magnetic field. On the other side (left) 
of the XDL, several unusually bright X-ray sources have also 
been detected. They are often referred to as /3 Ceti-like stars, 
based on the prototype (3 Ceti (Hiinsch et al. 1996, Maggio 
et al. 1998), which exhibits an extremely high X-ray lumi- 
nosity and a slow rotation rate (v sin i ~ 3 kms" 1 ). These 
observations are puzzling since it is difficult to reconcile the 
high X-ray emission with the slow rotation rate. Below we 
investigate the physical effects related to the accretion of a 
planet/brown dwarf in terms of the generation of magnetic 
field and X-ray emission. 

Suppose that an initially slowly rotating star undergoes 
an accretion event. As the planet spirals in, it will deposit an- 
gular momentum in the envelope and finally dissipate close 
to the core. However, the resultant profile of angular velocity 
ui(r) may not be uniform throughout the envelope but rather 

t A compact companion such as a white dwarf can be a relatively 
strong X-ray source if it accretes matter from the giant through an 
accretion disk. But such a companion would have been detected 
in their observations. 
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Figure 8. The parameters 7ce as a function of the radius for 
the three initial models (see text). In the expression for r^ ecsy , 
we took Mbd = i0Mj U p. The solid, dotted and dashed lines cor- 
respond to the initial models #A, #B and #C, respectively. 



may be steeper close to the dissipation region, where the 
gravitational potential increases abruptly. To quantify this, 
at least approximatively, we follow Livio and Soker (1988) 
and estimate the parameter 7ce defined by 



7CE 



T dc 



(4) 



where Td OC ay represents the decay timescale of the planet's 



orbit, and r m 



is the spin-up timescale of the envelope. 



The decay timescale is given by 
gM bd M(a) 



^dccay 



da/di 



1 



2a 



, V 



(5) 



where M(a) the mass interior to a and -Fd rag is the gravita- 
tional drag force that causes the spiralling-in. Its expression 
is approximately given by 

(6) 



drag 



where £ is a factor of order 2 — 5 for supersonic flows (e.g. 
Shima et al. 1985), R a ~ 2QM/V 2 is the accretion radius, 
V the relative velocity between the planet and the enve- 
lope and p the local density. Finally, the spin-up timescale 
7"s P in-up can be approximated by 

l(a)V 



a 2 F d 



(7) 



where 1(a) is the moment of inertia of the envelope interior 
to radius a. In the computations for our presented models, 
we assumed that the relative velocity V is Keplerian and 
impose £ = 4 as in Livio & Soker (1988). Figure |§| depicts 
the profiles of 7ce for the 3 initial models. Although these 
profiles are crude estimates, they suggest that a strong dif- 
ferential rotation is likely to take place at the base of the 
convective envelope, in the regions where 7ce < 1- The shear 
can then sustain a dynamo action and lead to a high X-ray 
emission (e.g. Regos & Tout 1995). Furthermore, the profiles 
of 7ce indicate that for less evolved configurations the outer 
layers are not spun-up, and therefore that it is possible to 
obtain X-ray emission from apparently slowly rotating stars 
as observed among /3 Ceti-like stars. It is difficult to estimate 
for how long the slow rotation will be maintained at the sur- 
face of the star since such a study would require to take into 



account the angular momentum transport in the convective 
envelope coupled with the presence of a strong magnetic 
field. Additionally the outcome depends on the evolution- 
ary status of the star. We expect however that the spin up 
the envelope will take longer in younger and less centrally 
condensed stars for which 7ce > 1 in a larger fraction of the 
envelope. In some configurations magnetic braking may not 
allow a significant spin-up of the outer layers. Some indi- 
cations that differential rotation may persist in the interior 
of stars are suggested by obs ervations. Indeed, in the HR 
diagram, the rotational (£4.2) and X-ray dividing lines are 
shifted with respect to one another. Thus, stars located be- 
tween these two lines present both a magnetic activity and 
a slow rotation rate. If this shift between the dividing lines 
corresponds effectively to an evolutionary sequence, this in- 
dicates either that the rotational braking timescale is shorter 
than the decay timescale associated with the dynamo activ- 
ity or, alternatively, that differential rotation is still main- 
tained in the stellar interior and sustains a dynamo activity. 

To summarize this section, the deposition of angular 
momentum by a spiralling-in planet/brown dwarf can pro- 
duce a strong differential rotation at the base of the convec- 
tive envelope. We then suspect that in this region a magnetic 
field can be generated via an a — Q dynamo and that this 
may ultimately result in an enhanced X-ray activity. Note 
that the magnetic field can also be used to brake the star. 
Finally, we would like to mention that the shear may trig- 
ger some insta bilit ies that could give rise to extra mixing 
processes (see 



4.4 The ejection of shells 

The increase in the mass loss rate during the accretion 
phase can lead to the formation of a shell around the star, 
and that shell can be expected to emit in the infrared 
(IR). More quantitatively, the mass loss rate can be in- 
creased by a factor of 2 to ~ 100 and a mass of the or- 
der 10- 6 M Q £ M cjcct £ 10- 3 M o can be ejected. These 
numbers depend however on the evolutionary status of the 
star, the duration of the accretion event and the accretion 
rate. Figure ^| demonstrates the very strong dependence of 
the mass loss rate on the accretion rate and on the initial 
model. The higher the accretion rate, the larger is the in- 
crease in the mass loss rate. Similarly, the more evolved the 
star, i.e. the lower its surface gravity, the easier it is to re- 
move material from its surface. Figure ^| also indicates that 
the rate of increase in Mi oss is faster for older stars. This is 
an illustration of the thermal response of the star since Mi oss 
is inversely proportional to the Kelvin Helmholtz timescale 
(A/i oss oc M/tkh)- Finally, we would like to emphasize that 
in the high accretion rate cases (cases B2 and C2), even a 
relatively small planet of a few Jupiter masses can increase 
the mass loss rate by about 2 orders of magnitude in a few 
hundred years! 

It is worth noting that in the estimates of the mass loss 
rate, the effects of rotation have been neglected. The depo- 
sition of the orbital angular momentum into the envelope 
can lower the effective surface gravity and facilitate mass 
ejection; it can also influence of geometry of the ejection. 
In particular, if the rotational velocity of the star exceeds 
some fraction of the break-up speed, then the ejected ma- 
terial could form an equatorially compressed outflow (e.g. 
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Figure 9. The evolution of the mass loss rate is shown for our 
different models. The solid lines refer to the standard evolution. 
In the lower panel, the dotted, short-dashed, long-dashed and 
dottcd-dashed lines correspond to cases AO, Al, A2 and A3, re- 
spectively. The same lines refer to cases BO, Bf, B2 and B3 in 
the middle panel and to cases CO, CI, C2 and C4 in the upper 
panel, respectively. The peaks in Mi oss in the upper panel are due 
to thermal pulses. 



Bjorkman & Cassinelli 1993, Owocki et al. 1994, Livio 1997). 
The peculiar circumstellar environment of U Equulei could 
be consistent with such a scenario (Barnbaum et al. 1996). 

The IR emission associated with giant stars is not easy 
to explain (Zuckerman et al. 1995), while in main sequence 
stars (Vega-like stars) the emission most likely arises from 
the dust associated with remnants of the star formation pro- 
cess or with small bodies like comets or planetesimals (e.g. 
Aumann et al. 1984, Telesco & Knacke 1991, Backman & 
Paresce 1993). In an extensive survey of luminosity class III 
stars, Plets et al. (1997) report that ~ 8% of the giants in 
their sample have IR excesses. They conclude that the IR 
emission is unlikely to be due to a present-day mass loss 
episode because, as they argue, (i) in a standard evolution, 
mass loss is predominantly efficient during a He shell flash 
so one would expect the stars to be more clustered at higher 
luminosities in the HR diagram and (ii) the dust emission 
is too cold for a recent ejection. They thus attribute the 
IR excess to a reheated dusty disk, as in Vega-like stars. 
Instead, we propose that the ejection of a shell resulting 



from an accretion event can account for the IR emission. 
Indeed, this scenario does not imply a clustering of stars at 
the tip of the RGB since the accretion event can, in prin- 
ciple, take place any time once the star leaves the main se- 
quence. Early engulfing could also account for the cold dust 
emission. Moreover, the IR excess is often too large to be at- 
tributed to a Vega-like phenomenon, since during the main 
sequence phase, most of the circumstellar debris will prob- 
ably be evaporated (Zuckerman et al. 1995). An important 
test of our scenario can be provided by measuring the 7 Li 
abundance for example, since this element should, at some 
level, accompany the accretion process (see §4.6). 
The detection of a detached shell around the J-type car- 
bon star Y Canum Venaticorum (Izumiura et al. 1996) is 
also interesting since this star is probably in an early evo- 
lutionary stage, on the RGB or on the He core-burning 
stage. The mass loss rate at the time of shell formation 
and the mass of the shell are estimated to be in the range 
(7-2O)lO" 6 M yr- 1 and (4- 1O)1O" 2 M , respectively. The 
present mass loss rate estimated from CO observations (e.g. 
Olofsson et al. 1993) is about 10 — 7 Mqvt -1 , two orders of 
magnitude smaller than the past value! Therefore, if the 
RGB status of Y CVn is confirmed, the existence of a mas- 
sive shell around this star may be a signature of an accretion 
event, since no other straightforward model can account for 
these observations. In this context, the massive shell and 
high mass loss rate at the time of ejection would suggest a 
very high accretion rate (M acc 10 -4 Mq yr -1 ). 

Shells have also been observed around many AGB stars 
(e.g. Olofsson et al. 1992, 1996, Bujarrabal & Cernicharo 
1994), but such structures can result from episodic high 
mass-loss events triggered by helium shell flashes (e.g. Olof- 
sson et al. 1990, Izumiura et al. 1997). It is therefore diffi- 
cult to identify the signature of an accretion event, although 
the presence of a companion is frequently invoked for the 
shaping of planetary nebulae (e.g. Soker 1997, Livio 1997). 
Finally, we should mention that in model CO the star expe- 
riences a series of thermal pulses before reaching the hori- 
zontal branch. From, the simulations we can estimate that 
during every interpulse, a mass of the order of 4 x 10" 4 M Q 
is ejected from the surface of the giant. These episodic ejec- 
tions would generate distinct shells, separated from each 
other by a few thousand years. However, the probability 
that a massive planet would be engulfed at the short-lived 
tip of the giant branch is rather small. 



4.5 Influence on the morphology of globular 
clusters 

The distributions of stars along the horizontal branch (HB) 
of globular clusters varies from one cluster to another. In 
some clusters, the population of HB stars exhibits an ex- 
tended blue tail while others contain only a red HB. It was 
first recognized by Sandage and Whitney (1967) that the 
metallicity is the first parameter influencing the morphology 
of the HB. Generally, more metal-rich clusters have red HBs, 
while metal-poor clusters have bluer HBs. But, for globular 
clusters of the same metallicity, the HB morphology can vary 
widely. Thus some "second parameter" other than metallic- 
ity must be affecting the HB morphology. Among potential 
second parameters, it has been shown that age (e.g. Lee et 
al. 1994), mass loss (e.g. D'Cruz et al. 1996), He mixing 
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(Sweigart 1997), CNO enhancement (Dorman et al. 1991), 
core rotation, or structural properties of the globular clus- 
ters such as concentration and density can also affect the 
HB morphology. Nevertheless the discovery of gaps in the 
blue tail of the HB (e.g. Walker 1992, Catelan et al. 1997, 
Ferraro et al. 1997) is still considered puzzling. 

Recently, Soker (1998b) suggested that planetary sys- 
tems can influence the HB morphologies. He proposed that 
the interaction of a RGB star with a giant planet/brown 
dwarf can affect the mass loss rate and thus the subse- 
quent location of the star on the HB. By depositing an- 
gular momentum and energy into the envelope, additional 
mass can be removed from the star during the RGB. The 
"peeled off" star will then have a higher effective tempera- 
ture and will appear bluer in the HR diagram (e.g. Castellani 
& Castellani 1993, Dorman et al. 1993). In particular, Soker 
(1998b) argued that, depending whether the planet survives 
the spiralling-in, collides with the stellar core or evaporates 
in the envelope, one can account for the 3 gaps observed by 
Sosin et al. (1997) in the blue tail of NGC 2808. 

Our study introduces a quantitative element into the 
scenario proposed by Soker. Our computations indeed indi- 
cate that a significant increase in the mass loss rate accom- 
panies the accretion process, especially if the star is evolved 
and the accretion rate is high. The simulations also show 
that, when the star has relaxed, the mass loss rate remains 
higher than in a standard evolution. We estimated that when 
a massive planet is accreted (not a brown dwarf), the mass 
loss rate is enhanced by a factor of 1.3 — 2 after the termi- 
nation of the accretion process. This is equivalent to setting 
Reimer's parameter tjr to a value close to 0.7 — 1.0 (in- 
stead of 0.5) in a standard evolution. If we now refer to 
the computations of Castellani & Castellani (1993) for such 
values for r/R, we conclude that stars accreting a planet at 
the bottom of the RGB, will actually be extremely blue on 
the horizontal branch. On the other hand, if the accretion 
occurs later, the effects are difficult to quantify because of 
the strong sensitivity of our computations to the accretion 
rate and the planet mass. The effects can be moderate, as 
in cases C2 or B2, or dramatic, as in cases B3 or C4. The 
possibility that a planet accretion would modify the evolu- 
tionary age of the star and thus affect the morphology of the 
HB is improbable. Indeed our initial models are relatively 
old (to > 12 Myr) and the duration of the red giant phase is 
relatively short compared to its age. Also, the accretion of 
a planet does not modify substantially the mass of the star 
therefore we don't expect the stellar evolutionary timescales 
to be significantly modified. 

Sweigart (1997) suggested that an extra mixing driven 
by internal rotation could transport helium from the HBS 
into the convective envelope. This would increase the RGB 
tip luminosity and hence the mass loss, leading to the for- 
mation of blue HB stars. Sweigart and Catelan (1998) also 
showed that only non-canonical mechanisms invoking deep 
helium mixing, internal rotation or a high initial helium 
abundance can reproduce the distribution of the blue HB 
stars in NGC 6388 and NGC 6441. However, these authors 
did not specify the source of the angular momentum required 
for the rotation scenario nor how the deep He-mixing pro- 
cess is activated. The accretion scenario can provide some 
support to these non-canonical processes, since it both offers 
a natural source of angular momentum for the core and it 



can be responsible for the generation of st rong instabilities 
as the planet collides with the core (§ 4.3 ) . Note also that 



some rapid rotators are observed in the blue side of the HB 
(e.g. Harris & McClure 1985). 

Concerning the gaps, recent observations (Ferraro et al. 
1998) indicate that they are not randomly distributed along 
the HB but rather located at similar positions in the color- 
magnitude diagram. Aside from the star-planet interaction 
hypothesis (Soker 1998b), it has been suggested that some 
poorly specified physical processes must operate during the 
RGB to increase the mass loss rate (Ferraro et al. 1998) 
or that gaps are a consequence of a bimodal or multimodal 
mass distribution of the globular clusters (e.g. Catelan et 
al. 1998). Let us therefore return to the planet accretion 
scenario and analyze the implications of this assumption. 
We have seen that depending on the evolutionary status of 
the star at the time of engulfing (i.e. on the initial planet 
separation), different amounts of mass can be lost and thus 
different positions in the HB can be reached. Assuming for 
the moment that planets form preferentially around solar- 
type stars (as it has been observed so far), the similarity in 
the positions of the gaps would suggest that planets orbit 
at preferred distances to the star. A theory suggesting that 
this is indeed the case exists (Notalle 1996, Notalle et al. 
1997 and references therein). If this is true, planets would 
be swallowed only when the solar-type star reaches a certain 
radius (or equivalently an evolutionary status), determined 
approximatively by the orbital separation of the planet. In 
this hypothesis, the accretion process would modify selec- 
tively the mass loss rate and thus the location of the star on 
the HB. 

However, the question of why the star-planet interaction 
process acts only sporadically and does not appear as a gen- 
eral phenomenon affecting the morphology of all the clusters 
would remain puzzling. To attempt to answer this question, 
one needs a better understanding of the processes involved 
in planetary formation. For example, metallicity and other 
globular cluster properties such as density may be impor- 
tant factors that can influence the formation of planets and 
they are still poorly understood. To conclude this section, 
the accretion of planets by giant stars can potentially influ- 
ence the mass loss rate and thus modify the position of the 
star on the HB. However, it is not clear at present that the 
gaps observed in the blue tail of some globular clusters are 
really the signature of such accretion events. 

4.6 The lithium rich G and K giants 

The detection of lithium in the spectra of G and K giants 
(e.g. Wallerstein & Sneden 1982, Brown et al. 1989, Grat- 
ton & D'Antona 1989, Pilachowski et al. 1990, Castilho et al. 
1995, de la Reza & da Silva 1995, Castilho et al. 1998) poses 
a problem concerning the origin of this element in these stars 
since, in standard evolution, lithium is entirely depleted dur- 
ing the first dredge up. Several mechanisms have been sug- 
gested to explain the Li enrichment and we refer to Gratton 
and D'Antona (1989) for a more complete discussion. Briefly, 
it has been proposed for example that Li could be preserved 
in the atmosphere of giant stars, but the low 12 C/ 13 C ra- 
tio observed in these stars (e.g. Pilachowski et al. 1990, de 
la Reza et al. 1995) indicates that deep mixing (and thus 
Li burning) has occurred. This represents therefore a failure 
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of this scenario. Gratton and D'Antona (1989) considered 
the possibility that the giant's envelope could be contam- 
inated by a nova exploding companion. However, the non 
detection of a white dwarf in the vicinity of Li-rich giants 
(e.g. de la Reza & da Silva 1995) weakens this hypothesis. 
The possibility that Li could be produced by the Cameron- 
Fowler mechanism (1971) has also been investigated. This 
mechanism, which consists of the rapid dredge up of 7 Be and 
Li production through the 7 Be(e~, v) 7 \A reaction, is known 
to operate in AGB stars of 4-6M© for which the convective 
timescale is short compared to the 7 Be destruction timescale 
(Sackmann & Boothroyd 1992). Unfortunately, this process 
is inefficient in low luminosity stars, unless one invokes some 
extra mixing processes (we will return to this point later). 
Finally, and we will develop this point below, several au- 
thors (e.g. Alexander 1967, Brown et al. 1989, Gratton and 
D'Antona 1989) proposed that the engulfing of a planet or a 
brown dwarf could be responsible for the high Li abundance. 

Based on IRAS observations of G and K giants, de la 
Reza et al. (1996, 1997) concluded that almost all the Li 
rich G and K giants have an IR excess compatible with the 
presence of a circumstellar shell. To account for these obser- 
vations, they constructed a model in which Li is produced 
inside the star by an extra mixing process, the "cool bottom 
burning" (Wasserburg et al. 1995). This process allows the 
fast transport of 7 Be from the HBS into the convective en- 
velope where 7 Be then decays to produce 7 Li. Furthermore, 
they assumed that the mechanism for Li enrichment is ac- 
companied by a sudden increase in the mass loss rate and 
thereby the ejection of a shell. Following the evolution of 
a thin expanding circumstellar shell, they could reasonably 
account for the IR emission and Li abundance distribution 
of G and K giants in a color-color diagram. However, de la 
Reza et al. remained relatively vague concerning the mecha- 
nism that is responsible for mass loss; they assumed that it is 
a consequence of the Li enrichment process. Instead, we pro- 
pose here that the accretion of a planet or a brown dwarf 
can naturally and consistently explain all of the observed 
features. We describe our scenario in the next paragraph 
and illustrate it in Fig. ^ which is an adaptation of Fig. 2 
of de la Reza et al. (1996). 

Figure [l(j represents a color-color diagram based on the 
IRAS density fluxes at 12, 25 and 60 /im. In this diagram, 
de la Reza et al. (1996, hereafter DDD) identified 3 regions, 
populated respectively by normal K giants with no IR ex- 
cess (region 7), Li rich giants with a large IR excess (region 
77), and relatively Li rich giants with a moderate IR excess 
(region 777). When the accretion process starts, Li is pro- 
gressively deposited in the stellar envelope and the mass loss 
rate increases. Rapidly, a shell is ejected and the giant moves 
quickly from region I to region II. In their shell model, DDD 
found that several hundred years are required to reach re- 
gion 77. This timescale gives an estimate of the duration of 
the accretion process. Then, Li starts to be depleted in the 
envelope while the shell expands and cools down. The 25/im 
flux decreases and the star now settles in region 777. DDD 
estimated that it takes approximatively 10 4 yr to move from 
region 77 to 777. Finally, in about 8 10 4 yr, the circumstellar 
shell dissipates, Li is entirely depleted and the star returns 
to region 7 as a normal K giant. Assuming a wind velocity of 
2 kms _1 and a mass loss rate of 5 10 -8 Mq yr _1 , DDD found 
that ~ 1O _5 M0 was ejected in the shell, in agreement with 
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Figure 10. The schematic evolution in an IRAS color-color dia- 
gram of a K giant star engulfing a planet or a brown dwarf. (By 
definition, [Ai — A2] = \og\iS\ 1 — logAiS\ 2 , where S\ is the 
IRAS density flux.) This figure was adapted from de la Reza et 
al. (1996) and is explained in the text. 



our models. In a second publication, de la Reza et al. (1997) 
completed their study by presenting new observations. They 
showed that Li-rich and Li-poor stars are mixed in the color- 
color diagram. They interpreted these observations as being 
due to differences in the Li destruction timescales. In our 
model, this is explained naturally if the amount of Li de- 
posited into the star, or equivalently the mass of the planet, 
is different from one star to another. 

The mechanism that leads to the Li enrichment must 
also account for its rapid depletion, which takes place in 
~ 10 5 yr- In our computations without diffusion, Li is main- 
tained at the surface of the star but Charbonnel (1995) 
showed that this element can be rapidly depleted in red gi- 
ants due to non-standard mixing processes. In particular, 
the model developed by Zahn (1992) showed that the ef- 
ficiency of the mixing is mainly determined by the loss of 
angular momentum via a stellar wind and, even in the ab- 
sence of a high mass loss rate, the process takes place where 
there is a strong gradient in ui. He also pointed out that 
additional mixing is expected in shell burning regions. Since 
all of these conditions are satisfied in the context of the ac- 
cretion process (increase in the mass loss rate, deposition of 
angular momentum and a str ong gradient of u) in the dissi- 
pation region as described in £4.3) and as a part of the fact 
that the star is on the RGB (the presence of the HBS), we 
have every reason to believe that the mechanism proposed 
by Zahn (1992), and that successfully explains some pecu- 
liar anomalies in red giant stars (Charbonnel 1995), could 
be responsible for the fast Li depletion following a short ac- 
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cretion phase. Moreover, observations of the 12 C/ 13 C ratio 
in some Li rich giants (da Silva & de la Reza 1995) indicate 
that some extra mixing is operating since the carbon isotope 
ratio in these stars is lower than the standard value. 
Our scenario does encounter some difficulties in explaining 
the very high Li abundances [loge(Li)<; 2.8] found in some 
giants (e.g. Gratton & D'Antona 1989, de la Reza & da Silva 
1995, Fekel et al. 1996). This would require the deposition 
of a huge amount of lithium that can only be achieved either 
by the accretion of numerous planets such that M acc » M elw 
or if the Li abundance in the brown dwarf is extremely high 
(X acc £ 10 x X cnv for Af acc ~ 0.1M Q ). The accretion sce- 
nario does not provide a straightforward explanation but we 
can imagine that when the planet reaches the central region, 
the dissipation process may generate some instabilities (due 
to the shear) that could, in turn, trigger the "cool bottom 
burning" invoked by de la Reza et al. (1996, 1997). The ac- 
cretion process would then lead to some Li production (as 
in more massive AGB stars). 

De Medeiros et al. (1996a) reported that the majority 
of Li K giants are normal rotators (usini ~ 2 — 3 kms -1 ) 
although some exceptions exist [e.g. HD 233517 (Fekel et 
al. 1996), HD 9746, HD 25893, HD 31993, HD 33798, HD 
203251 (Fekel & Balachandran 1993), IE 1751+7046 (Am- 
bruster et al. 1997)]. However, as pointed out in §4.2, high 
angular velocities (~ 10% percent of the Keplerian veloc- 
ity) can only be achieved if a sufficiently massive planet 
(Mhd > 5Mj U p) is accreted. For a typical 1.0M© RGB star, 
with a radius R ~ 20_Rq , the deposition of the orbital angu- 
lar momentum of a 5 Mj up planet will "only" increase the 
star's rotational velocity to ~ 5.5 kms -1 [Eq. (0)]. Conse- 
quently, as long as the mass of the engulfed body is low (a 
planet rather than a brown dwarf), the observed angular mo- 
mentum of Li rich giants is compatible with our scenario. We 
should also note that due to the increased mass loss, mag- 
netic braking could also be operating (e.g. Leonard & Livio 
1995). Finally, let us point out that all the Li rich giants 
are single stars, hence according to the prevailing common 
wisdom, they could possess planetary systems. 

To summarize this section, we propose that the accre- 
tion of a planet or a brown dwarf is responsible for both the 
IR excess and the high Li abundance present in almost all 
single Li rich G and K giants. 



5 DISCUSSION AND CONCLUSION 

In this paper, we investigated the effects of the accretion of 
a planet or a brown dwarf by a solar type star located on 
the red giant branch. Our computations show that the de- 
position of mass (potential energy) produces a substantial 
expansion of the envelope and, for higher accretion rates, 
nuclear burning at the base of the convective envelope can 
take place. The simulations show that the results are a sen- 
sitive function of the accretion rate and the structure of 
the star. Generally, the higher the accretion rate and/or the 
more evolved the star, the larger is the expansion. After the 
end of the accretion process, the star relaxes and resumes 
a standard evolution. However, in one simulation (case CO), 
the star avoided the He flash but instead developed thermal 
instabilities. For this model, after a phase of thermal pulses, 
He burning was finally ignited at the center of the star. 



We discussed the possible observational signatures 
of the accretion of a planet /brown dwarf. In particular 
(1) we showed that the rotational velocity of the star 
can be substantially increased, provided that the accreted 
planet/brown dwarf is sufficiently massive. In this respect 
the accretion of a brown dwarf could account for the fast ro- 
tation of FK Comae stars and more generally, the deposition 
of angular momentum via this process can explain unusually 
fast rotation rates. (2) Investigating the X-ray properties of 
giant stars, we showed that a strong differential rotation may 
take place at the bottom of the convective envelope. The 
presence of shear in this region can turn-on the dynamo ac- 
tivity and this may ultimately generate X-ray emission. We 
showed that under some conditions the surface layers of the 
star could rotate slowly while a rapid rotation is maintained 
in the core. (3) We also showed that the accretion scenario 
can lead to the ejection of shells whose mass is in the range 
10- 6 M Q <, M cjcct <, 1O" 3 M . The IR excesses observed in 
~8% of luminosity class III stars could in part be due to this 
process. (4) We argued that the increase in the mass loss rate 
can also have indirect effects on the morphology of the hor- 
izontal branch and in this respect, the accretion process is a 
potential second parameter candidate. (5) We showed that 
the deposition of planet material in the convective envelope 
can efficiently increase the 7 Li surface abundance and more 
generally enhance the stellar metallicity of giant stars. (6) 
Finally, on a related topic, we proposed that the accretion 
of a planet/brown dwarf can explain consistently both the 
IR excess and the high Li abundance observed only in single 
red giant stars. 

In terms of physical processes, we suggested that 
planet/brown dwarf accretion may provide the triggering 
for a variety of processes. For example, the disruption of 
the planet and the associated instabilities could increase 
the He abundance in the giant's envelope and thus affect 
the morphology of the HB as proposed by Sweigart (1997, 
§1.5). They could also be responsible for the production of 
lithium th roug h the mechanism proposed by Wasserburg et 
al. (1995, §4.6). Finally, the rapid contraction of the star at 
the end of the accretion process could generate some pulsa- 
tions. 

Observations indicate that several of potential different 
signatures of an accretion event, e.g. a fast rotation rate, a 
significant X-ray emission, a high Li abundance or a large 
IR excess, do not always occur simultaneously. As we have 
shown, this is not surprising. First, depending on the mass 
of the planet, the outcomes could be different. For exam- 
ple, a small planet that evaporates in the envelope before 
reaching its base will perturb the structure only slightly. In 
such a case, no shell will be ejected but the deposition of 
angular momentum and matter in the surface layers could 
still produce a relatively rapid rotation and a high lithium 
abundance, as observed in some chromospherically active gi- 
ants (Fekel & Balachandran 1993). Secondly, the timescales 
associated with the various phenomena (IR excess, X-ray 
activity, rotational braking, particle transport) can be dif- 
ferent, and thus, not all the signatures will be observed at 
the same time. For example, if the spin-down timescale is 
longer than the timescale associated with the dissipation of 
the ejected shell, one can account for the absence of IR ex- 
cess in some fast rotators. Thirdly, the consequences of the 
accretion process depend on the evolutionary status of the 
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star. Giant stars located at the bottom of the RGB have 
much deeper convective envelopes than near the tip of the 
giant branch. Thus, if we consider for example the deposi- 
tion of Li in the envelope, it is more difficult to obtain high 
Li abundances in younger red giants (it is also be more diffi- 
cult to eject a shell). Finally, our simulations show that the 
results are sensitive to the physics of the dissipation process, 
and more specifically, to the accretion rate and the location 
of the accretion process. At present, both these physical pa- 
rameters are rather uncertain, and a better understanding 
of the dissipation process is definitively required. 

As a last point, we can attempt to use the expected 
observational signatures to infer some estimate for the fre- 
quency of planets around solar-type stars. In a search for 
lithium in 644 G and K giants, Brown et al. (1989) found 
that 2% in the stars of their sample have loge(Li)> 1.5 and 
4% have loge(Li)> 1.3. If we adopt for our definition of Li 
rich giants the value loge(Li)> 1.2, we find that ~ 8% of 
the stars surveyed by Brown et al. belong to this group. In a 
more recent survey, Castilho et al. (1998) report the detec- 
tion of 5 Li rich and 6 moderately Li rich giants among 164 
stars, indicating a proportion of Li rich giants of ~ 3 — 7%. 
The study of Plets et al. (1997) indicates that 8% of lumi- 
nosity class III giants have unexpected IR excess. Finally, 
based on rotational velocity measurements of ~ 900 G and 
K giants, De medeiros et al. (1996a) found that <y 5% of 
late G and K giants located beyond the rotational dividing 
line are fast rotators with v sin i > 5 kms^ 1 . Thus, if the 
IR excess and/or high Li abundance and/or unusually high 
rotational velocity among G and K giants result from the ac- 
cretion of a planet/brown dwarf, the above numbers suggest 
that at least 4-8% of the single, solar-type stars must have a 
low mass companion, either a brown dwarf or a planet. If we 
were to take seriously the frequency of planets inferred by 
Soker (1997) to account for the shaping of planetary nebulae 
(~ 55% of all progenitors of masses <, 5M Q ), our deduced 
frequency would indicate that either a large fraction of the 
planets do not interact during the red giant phase or, more 
probably, that the observational signatures of the accretion 
event are very short-lived compared to stellar evolutionary 
timescales. Indeed, because of the large radii reached by so- 
lar type stars during the RGB, some planets will most likely 
be engulfed during the early giant phase rather than on the 
AGB (although Jupiter will not be engulfed by our own sun). 
On the other hand, if we compare the timescales associated 
with both the Li and IR excesses (i; 10 5 yr) to the red giant 
lifetime (between 6 10 6 and 5 10 8 yr for giants with masses 
between 1.0 and 2.5M Q ), we conclude that the probability 
of detecting an accretion event is such that we are likely to 
miss most candidates. 
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